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The scalar relativistic version of an 
accurate first principles full potential 
self- consistent linearized muffin tin 
orbital (LMTO) method has been employed for 
describing the physical properties of the 
parent system of the high_Tc oxide 
superconductors, i.e., CaCu02. The presently 
employed modified version of the LMTO method 
is quite fast and goes beyond the usual 
LMT0_ASA method in the sense that it permits 
a completely general shape of the potential 
and the charge density. Also, in contrast to 
LMT0_ASA, the present method is also capable 
of treating distorted lattice structures 
accurately. The calculated values of the 
lattice parameters of pure CaCu02 lie within 
3% of the experimentally measured values for 
the Sr-doped system Ca( .86)Sr( .14)CuO<2) . The 
computed electronic structures and the 
density of states is quite similar to those 
of the other oxide superconductors, except of 
their three- dimensional character because of 
the presence of strong coupling between the 
closely spaced Cu02 layers. The van Hove 
singularity peak appears slightly below the 
Fermi level and a small concentration of 
oxygenation /or/ substitutional doping may 
pin it at the Fermi level. The calculated 
frequencies for some symmetric frozen phonons 
for undoped CaCu02 are quite near to the 
measured data for the Sr-doped CaCu02. 


I. INTRODUCTION 

The present LMTO method is seen to produce 
the electronic structure, cohesive energy, lattice 
constants, elastic constants, phonon frequencies, 
mode Grtieneisen and strain parameters for the simple 
systems like Si, C etc [1J. Very recently, the method 
has been successfully applied also for the III-V and 
I I -VI semiconducting compounds like AlAs, CdTe, GaSb, 
ZnSe, ZnTe, ZnS etc [2,3]. The influence of 
structural relaxation of the atoms on the 
valence-band off-set at the lattice matched 
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interfaces of II-VI and III-V semiconductors 
ZnTe/GaSb( 110 ) and the lattice mismatched interface 
ZnS/ZnSe( 001) has been investigated [41. Very 
recently the van Hove singularity scenario and hole 
concentrations at different pressures for the 
Hg-based oxide superconductors have been investigated 
[ 5 ] . 


A compound CaCu02 is obtained if the number of 
Cu02 planes is allowed to increase indefinitely (n -> 
oo ). Thus, CaCu02 is the infinite number of Cu02 
layers limit of the above superconductors and is seen 
to be insulating 16} . Tc might increase with n. The n 
= qd material appears to be more difficult to dope. 
However, more recently (CaSr)Cu02 has been doped 

successfully and an electron type superconductivity 
has been seen [7]. In practice, the n = od structure 

is not formed with only Ca ion on A site, a material 

with some Ca substitution by Sr, 

Ca ( 0 . 86 ) Sr ( 0 . 14 )CuO( 2 ) has been prepared. The 
Sr-raixed compound is seen to be narrow band 
semiconductor. For details of the LMTO method we 
refer to the earlier papers [1-4]. 

II. CALCULATION AND RESULTS 

The CaCu02 structure as shown in Fig.l has a 
space group P4/mmm (D4hli{). The atomic positions in 
the unit cell are Cu: (0,0,0); 

0(1) : ( .5a, 0,0) ; 0< 2 ) : ( 0 , . 5a , 0 ) ; Ca: ( .5a, .5a, .5c) . 

For the Sr-doped 
sample the experimental 
lattice parameters are 
a = 3.86 A and c = 3 . 2 A° . 
The bond lengths are 
Cu-0<1) = Cu-0 ( 2 ) 

1.93 A° and Ca-O(l) = 
Ca-0( 2 ) = 2.50 A°. 

In the present 
calculation, the basis 
employed for making 
expansions of the 
products of the LMTO 
envelops included 

functions with 154 
and of energies -0.01, 
-1.0 and -2.3 Rydberg 
and with decays given 
by X =-l and -3 
Rydberg. The set will 
include 50 functions 



FIG.l Unit cell for the 
crystal structure 
of CaCu 02 • 
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for each atomic sphere. The local density potential 
of Hedin and Lundquist has been employed. An absolute 
convergence to better than 1.5 mRy/atom is obtained 
with spd basis of 22 LMTO' s/atom for each atom 
including 0. The number of atoms in unit cell were 
taken to be four. 

The muffin tin (MT) spheres were chosen to be 
slightly smaller than touching and the radii for Cu, 
0 and Ca were taken as 1.88, 1.64 and 2.99 atomic 
units(a.u.) f respectively. The states Cu( 3d, 4s , 4p) , 
0(2s,2p) and Ca(3p,3d,4s) are taken as the valence 
band electrons. An empty sphere has been included at 
(0,0,. 5c) in the calculation. 



FIG. 2 Variation of crystal energy 
(mRy) with the ratio of unit 
cell volumes, Vo and V as 
described in the text. 



FIG. 3 Energy variation with the 
interplaner spacing. 


Variati- 
ons both in 
the lattice 
parameter 'a' 
and the ratio 
c/a were 

investigated. 
The results 
are shown in 
Figs- 2 & 3, 

respectively . 
One finds 

energy 

minimum for 
a volume of 
the unit cell 
V = 0.925 Vo 
and c/a = 
1.00, where 
Vo is the 
unit cell 

volume for a 
= 3.86 A (in 
fact, these 
lattice 
parameters 
correspond to 
the 

experimental 
values for 
the Sr doped 
CaCu02 ) . So 
the values of 
a and c are 
quite close 
to those for 
Sr- doped 

material. The 
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calculated value of ' c' is within 3% of the 
experimental value. The calculated value of bulk 
modulus is 1.99 Mbar. 

The charge transfer from the muffin tin spheres 
(i.e., the difference between the atomic number and 
charge lying within the MT sphere ) of Cu, 0 and Ca 
are 1.65, 0.56, and 0.62 electron charges, 

respectively. The total electron charge per unit cell 
outside the MT spheres is 3.39e. 

CA) Electronic Structure 

The calculated band structure along some 
symmetry directions is presented in Fig. 4. The 
various symmetry points of the Brillouin zone in 
units of 2rr/a are r( 0,0,0); X (.5,0,0); M(.5,.5,.6); 

Z(0,0,.6); R( .5,0, .6)and A(.5,.5,.6). As the spacing 
between the two successive Cu02 layers is smaller 
than the intraplaner Cu-0 seperation within the 

Cu02 , the 

bands along 
the 

c-direction 
show 

appreciable 
dispersion 
although 
smaller than 
the 

dispersion 
seen within 
the x-y plane 
(or Cu-02 
layer) . The 
Fermi level 
is set at 
origin of 

energy. 

FIG. 4 Dispersion curvs for the 

CaCu0 0 system. The lowest 

* 11 bands 

lying in the 

energy interval range of about 9 eV near and below 

Fermi level arise mainly from the Cu(3d) and 0(2p) 

states. This spread (=* 9 eV) of energy bands is 

similar to that calculated by Mattheiss and Hamann 
[81 and by Singh et al [9] using LAPW method and by 
Korotin an Aniscimov [10] using the LMTO-ASA method. 
The uppermost band crossing the 2 Fejmi level 
corresponds to the anti-bonding Cu-d[x -y ]:0-p[x,y] 
states. The lower parts of the Ca(3d) bands lying 
above Fermi level overlap the antibonding 
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FIG. 6 Projected DOS for Cu, O 
and Ca atoms with s, p & 
d orbitals contribution. 


near r and A 
symmetric points. 

In order to 
see the nature of 
the dispersion in 
another plane 

shifted in the c - 
direction, we have 
also included in 
Fig. 4 the 

dispersion curves 
in the 

2rr/a[0,0,0.6] 
direction. There 
appear saddle 

points at -0.81 eV 
at X point and at 
-0.19 eV at R 
point. 

In the 

calculation of the 
total density of 
states (DOS), a 
sampling method 
for 196 k points 
in the irreducible 
part of the 

Brillouin zone 

with a Gaussian 
energy broadening 
of 0.002 Ry was 
employed. The 

calculated DOS for 
the self 
consistent 
calculation is 
shown in Fig. 5. 
The DOS shows the 
main features 

which are quite 
similar to those 
of the other oxide 
high_Tc cuprates . 
The projected 

densities at 

various atoms have 
also been shown in 
Fig. 6. At Fermi 
level, the main 
contributions 
arise from Cu(3d) 
and 0(2p) states. 
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The high 
degree of the nesting 
of the Fermi surface is 
likely to give rise to 
singularities in the 
generalized k-dependent 
electronic 

susceptibility of the 
two-dimensional system 
and may result in the 
electronically-driven 
instabilities such as 
the incommensurate 

charge or spin density 
waves . 

B. Phonon Frequencies 

The variation of the 
internal energy of the 
solid with the various 
types of static 

deformations in the 

small displacement 

limit has been employed 
for calculating the 

frozen phonon 

frequencies at the 

r<k=0) point of the 

Brillouin zone. In 

CaCu02 in all, there 
are six optical modes 
out of which five are 
infrared active and the 
remaining sixth mode comprising of the out of (Cu02) 
plane vibrations of oxygen atoms in opposite 
direction is silent. The frequencies of 

these vibrations alongwith their eigenvectors are 

presented in Table I. The frequencies have 

been calculated by determining the harmonic force 
constant after fitting the energies of the distorted 
structure for several amplitudes with a polynomial 
containing terms up to third degree. In order to see 
the reliability of the estimates of these calculated 
values we have also included in the table the 
measured values for the Sr-doped CaCu02 . The 

frequencies for undoped sample are quite close to the 
observed values for the Sr-doped sample except for Eu 
mode which involves the motion of all the atoms. 


Chart* Denaity of CaCuO* 
in (Of— *1,1) plana 



FIG. 8 Valence charge dens- 
ity contors in the 
(O f -l f l) plane in 
interval of 0.015 
e/ (a.u. ) up to a 
maximum value of 

0.12 e/(a.u.) . 
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CaCu0 2 


Fermi Surface 

r x r 



FIG. 9 Cross-sections of the Fermi surface in the 
extended Brillouim zone scheme in two planes 
(a) kz = 0.0 and (b) kz = 0.603 for CaCu02. 


III. CONCLUSIONS : 

The LMTO method which is comparatively 
simple and faster than the LAPW method is able to 
predict one electron energy spectrum of CaCu02 quite 
well. Values for the bulk modulus and the frozen 
phonon frequencies have been otained. The phonon 
frequencies for the undoped CaCu02 are quite near to 
those measured for the Sr-doped sample 


561 




TABLE I. Frequencies of the various phonon modes 
for CaCu02 in cm - 1 



ACKNOWLEDGEMENTS 

The financial assistance from University Grants 
Commission, New Delhi and Department of Science and 
Technology, New Delhi are acknowledged. 


562 













REFERENCES 


1. M. Methfessel, Phys. Rev., B38 , 1557 (1988); M. 

Methfeasel C.O. Rodriguez and O.K. Anderson Phys. 
Rev., B40, 2009 (1989). 

2. Bal K. Agrawal and Savitri Agrawal, Phys. Rev., 
B45, 8321 (1992). 


3. Bal K. Agrawal, Savitri Agrawal and P.S. Yadav, 
preprent (1994). 

4. M. Methfessel, Bal K. Agrawal & M. Scheffler, 
Proc. Int. Conf. of Physics of Semiconductors, 
1 990 , Greece; Bal K. Agrawal and Savitri Agrawal, 
preprent, (1994). 

5. Bal K. Agrawal and Savitri Agrawal, Physica C, 
(1994) in press; High Tc update, 1 May, (1994). 

6. T. Siegrist, S.M. Zahurk, D.W. Murphy, & R.S. Roth, 

Nature 334. 193 Cl 9883. 


7. X. Smith et al . , Nature 351, 549 (1991). 


8. L.F. Mattheiss and D.R. Haraann, Phys. Rev. B40 
2717 (1989) . 

9. D. Singh, W.E. Pickett and H. Krakacur, Physica C, 
162-164 , 1431 (1989). 

10. M . A . Korotin and V.I. Aniscimov, Mater. Lett. 10, 
28 (1990). 

11. G. Burns et al, Phys. Rev. B40 , 6717 (1989). 


563 



